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Abstract 

The calculated electron mobility for a graphene nanoribbon as a function of applied electric field 
has been found to have a large threshold field for entering a nonlinear transport regime. This field 
depends on the lattice temperature, electron density, impurity scattering strength, nanoribbon 
width and correlation length for the line-edge roughness. An enhanced electron mobility beyond 
this threshold has been observed, which is related to the initially-heated electrons in high energy 
states with a larger group velocity. However, this mobility enhancement quickly reaches a maximum 
due to the Fermi velocity in graphene and the dramatically increased phonon scattering. Super- 
linear and sub-linear temperature dependence of mobility seen in the linear and nonlinear transport 
regimes. By analyzing the calculated non-equilibrium electron distribution function, this difference 
is attributed separately to the results of sweeping electrons from the right Fermi edge to the left 
one through the elastic scattering and moving electrons from low-energy states to high-energy ones 
through field-induced electron heating. The threshold field is pushed up by a decreased correlation 
length in the high field regime, and is further accompanied by a reduced magnitude in the mobility 
enhancement. This implies an anomalous high-field increase of the line-edge roughness scattering 
with decreasing correlation length due to the occupation of high-energy states by field-induced 
electron heating. 
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I. INTRODUCTION 



The engineering achievement of isolating graphene sheets-"- from graphite has inspired 
many studies aimed at understanding basic underlying physics-i^ as well as finding possi- 
ble applications to carbon-based electronics-. The linear transport of charge carriers in a 
graphene layer, has received a considerable amount of attention. ^-^^ Recent reports on the 
successful fabrications of ultra-fast graphene transistors-i^ and photodetector have further 
pushe d this research frontier into the fields of electronics and optoelectronics. However, 
similar studies of linear transport in graphene nanoribbons (GNRs)has only been given 
relatively little attention so far.— ^-i^ 

Early theoretical studies— on electron transport in graphene nanoribbons were re- 
stricted to the low- field limit, where a linearized Boltzmann equation was solved within 
a relaxation-time approximation. In this paper, the non-equilibrium distribution of elec- 
trons is calculated exactly by solving the Boltzmann equation beyond the relaxation-time 
approximation for nonlinear electron transport in semiconducting graphene nanoribbons. 
Enhanced electron mobility from initially-heated electrons in high energy states is antici- 
pated. An anomalous increase in the line-edge roughness scattering for large electric fields is 
obtained with decreasing roughness correlation length due to the occupation of high-energy 
states by field-induced electron heating. Although the numerical results in this paper are 
given only for semiconducting graphene nanoribbons, the derived formalism here can also 
be applied to metallic graphene nanoribbons by formally taking the limit W ^ oo. The 
semi-classical Boltzmann transport equation is expected to be applicable to the diffusive 
band-transport regime with relative smooth edges for graphene nanoribbons, instead of the 
hopping and tunneling between localized states with rough edges. 

The rest of the paper is organized as follows. In Sec. [Tll we exactly solve the semi-classical 
Boltzmann equation for low-temperature electron transport in semiconducting graphene 
nanoribbons by including impurity, line-edge roughness and phonon scattering effects at a 
microscopic level. Based on the calculated non-equilibrium distribution as a function of wave 
number along the ribbon, we present detailed numerical results for electron mobility as a 
function of either the applied electric field or the lattice temperature for various electron 
densities, impurity scattering strengths, nanoribbon widths and correlation lengths for line- 
edge roughness. The conclusions deduced from these results are briefly summarized in Sec. 
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II. MODEL AND THEORY 

We restrict ourselves to single subband transport. Low electron densities, moderate tem- 
peratures, ionized impurities and line-edge roughness are assumed— i^^ii^. Therefore, neg- 
ligible effects of electron-electron—, optical phonon— , inter-valley and volume-distributed 
impurity scattering— will be ignored. Typical armchair nanoribbon (ANR) configuration is 
shown in Fig.[Tl The electron-like dispersion equation branch for ANR can be represented 
on a /c— mesh as- 

{kj , metallic 

/ • (1) 

y k'j + (ir /3W) , semiconducting 

Here, i'f = 10*cm/s is the Fermi velocity in graphene. The wave numbers kj = 
[j — {N + l)/2] 6k are given on the discrete mesh j = 1, 2, . . . , N for large odd integer 
A^. 6k = 2 k^a_^/{N — 1) is a small mesh spacing, /cmax is chosen so that scattering induced 
population of higher elctron-like branches may be neglected (see Fig.[T]). The minimum in 
the dispersion curve corresponds to the central mesh point j = M = [N + l)/2. Also, 
W = (A/" + 1) ao/2 is the width of the ribbon expressed in units of the size of the graphene 
unit cell = 2.6 A and number of carbon atoms across the ribbon A/". According to the 
dispersion relation in Eq. ([1]), the electron group velocity Vj for metallic nanoribbons on the 
first subband is determined by the Fermi velocity up, while for semiconducting ANRs, it 
assumes the form vj = vp {hupkj/ej). We shall also assume that the electron-like branch is 
filled np to \kj\ = kp at zero temperature with the Fermi wave number and energy given 
by kp = 7rni£)/2 and ep = e{kp), which can be applied to both metallic and semiconduct- 
ing graphene nanoribbons by choosing the zero energy at the subband edge. For a chosen 
temperature T and chemical potential /iq, the linear electron density in ANR follows from 

N 

niD = 6k/-n J2 [exp {{cj - fio)/kBT) + 
i=i 

The wave function "^jix, y) corresponding to the energy in Eq. ([T]) must satisfy hard-wall 
boundary conditions'^ ^j(O) v) = ^i(W^7 v) = 0- This can be assured by choosing the wave 
function as a mixture of those at K = (2n/3aQ, 27r/-\/3ao) aiid K' = (— 27r/3ao, 27r/-\/3cto) 
points as- 
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^A^^ y) = 



^i{2TT /3ao—K)x 



at K point 



(3) 



2LW 



-i{2TT /3ao—K)x 



at K' point 



Here, L is the quantization length of the ribbon. For semiconducting ANR k, = 7t/3W ^ 
27r/3ao is the quanta of the transverse wave vector, (pj = tan"-*^ [kj/ k) is the phase separation 
between the two graphene sublattices. For metalhc type of = ribbon, we have /t = and 
the phase assumes only ±7r/2 values. 

Given the wave function in Eq. ([2]) and neglecting effects due to inter-valley scattering- 
one can calculate the scattering from any given potential V{x^ y). The impurity and phonon 
induced inter-valley scattering has been ignored based on the assumption that the required 
momentum transfer is relative large in comparison with the low-temperature phonon ener- 
gies as well as the suppressed effective scattering cross section of both volume and surface 
impurities for very large value |K — K'|. Corresponding interaction matrix elements are 



w 



w 



Vi. 



dx I dy^*{x, y) V{x, y) ^^(a;, v) = ^ I dx 



dy {^:Vtlj,+tlji*V^'j) . (4) 



— oo ~oo 

Our model accounts for the net scattering potential made from the three components V = 

The longitudinal mode considered in this paper induces higher deformation potential 
than the out-of-plane flexural modes, and we neglect the transverse flexural modes here as 
a leading-order approximation. In this approximation, the inelastic scattering is provided 
only by longitudinal acoustic phonons. The perturbation potential induced by them is given 
byiS^ 



V^'iy) 



n 



2pLWuAL 



DalQv e 



(5) 



where the phonon distributions are n~ = [exp (huAL/^BT) — = 1 + n~, T is 

the temperature, ual = qy is the phonon frequency, Dal ~ 16 eV is the deformation 
potential, p ~ 7.6 x 10~^g/cm^ and ~ 2 x lO^cm/s are the mass density and sound 
velocity in graphene. Also, one has to assure momentum conservation during scattering, 
i.e., Qy = ki kj. 

Elastic scattering is due to the roughness of the ribbon edges as well as in-plane charged 
impurities. To model the first we shall assume that the width of the ribbon has the form 
W{y) = W + 6W{y), where the edge- roughness is assumed to satisfy Gaussian correlation 
function {6W{y) 6W{y + Ay)) = dlP' exp[— (Ay/Ao)^] with 56 ~ 5 A being the amplitude and 
Aq ~ 50 — 200 A being the correlation length of the roughness. The second mechanism is 
given by impurities located at (xq, 0, 0) and distributed with a sheet density n^u. Each of the 
point impurities provide the scattering Coulomb potential. Since the momentum difference 
between two valleys is rather large, |K — K'| = 47r/3a, only short-range impurities (such as 
topological defects) with a range smaller than the lattice constant causes inter- valley scatter- 
ing. In contrast, for wide enough ribbons with uniformly distributed long-range impurities 
(such as charged defects) the scattering processes is restricted to intra-valley scattering—"^. 
Here we neglect the short range scatterers in favor of the charged impurities thus making 
the scattering matrix diagonal with all elements given by Eq. (7). Such simplification has 
been invoked in order to keep the conventional part of the Boltzmann transport equation as 
simple as possible and focus on the nonlinear counterpart. Corresponding perturbations for 
roughness and impurity scattering are given by 



V-^-iy)^'-^.?..,. (6) 



e2 



V^"-P(a;, y) = , (7) 

where is the average dielectric constant of the host. These potentials give rise to the net 
elastic scattering rate h/rj on the /c— mesh 
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+ -LER ' (8) 




imp ^LER ' 



70 (a) [1 + c°<2</),)] 



(9) 



where 70 denotes the impurity scattering rate at the Fermi edges and 71 = 2 ( '^'^2^ )^ ^ is 
the scattering rate due to edge roughens. Here, we assume that the impurities he within 
the graphene layer. The scattering rate will be modified for charged impurities in a sub- 
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seen from the effective scattering cross section S{ky, W) de- 



strate, which can be easi 
fined in Eq. (28) of Ref. [16|. Although this effect can be formally taken into account 
by varying the value of 70 in Eq. ([9]), the leading scattering mechanism is still provided 
by in-plane impurities. All the scattering potentials are screened by free carriers. Since 
we are limited to a single subband, the screening is described by the dielectric function. 
The inelastic scattering is shielded by the static Thomas-Fermi dielectric function in its 
general form— e^Fd^j' ~ ^il)- The screening of elastic scattering potentials is given by 
eTF ~ 1 + (e^/n'^eoerhh'F) under the metalhc limit {2kFW ^ 1) with ~ 3.9. We employ 
the static screening to both impurity and phonon scattering with electrons for the system 
with a relative large damping rate and small Fermi energy. The static screening tends to 
overestimate the scattering effect.— The dynamic screening effect on the electron scatter- 
ing has been quantitatively evaluated within the RPA in a similar quantum-wire system,— 
and it becomes important whenever the damping energy is much smaller than the Fermi 
energy. The Fermi-golden rule has been employed to treated the leading contributions of 
both elastic and inelastic scattering. The multiple scattering effects are neglected due to low 
temperatures, low impurity densities and small line-edge roughness. Under application of a 

small electric field J-q along the GNR (see Fig.[T]) the stationary current / = (eSk/n) ^ fj vj 

j 

establishes itself. It can be characterized by the electron mobility fie = I/eniDTo, where we 
have introduced the carrier distribution function partitioned as f) = ff^ +gj ~ fj^^+bj r™^ 
with fj^^ = {1 + exp [{sj — fic) /^bT]}~^ being the equilibrium Fermi-Dirac function, /ic be- 
ing the chemical potential of electrons, and bj = evj {d ff^ / de j) J^q being the drag. Such 
approximation for the distribution function is known as relaxation time approximation and 
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encompasses all (elastic and inelastic) scattering mechanisms into r™^. The consequence 
of such an approximation is the electric field independent mobility and conductivity. How- 
ever, by increasing J-'q one may enter nonlinear regime when the conductivity and mobility 
become a function of the applied electric field. Studying of ANR in such nonlinear regime 
is the main subject of this work. 

Formally, the deviation from the equilibrium Fermi distribution under strong electric field 
is described by the set of nonlinear reduced Boltzmann transport equations— 

^=^.- E (10) 

In this notation, g'jit) = gj{t) — gM{t) is the reduced form of the dynamical non-equilibrium 
part ^® of the electron distribution function. The reduced form accounts for particle number 

N 

conservation condition, i.e., ^ gj{t) = 0. The nonlinear reduced Boltzmann transport 

equations have been applied to study the electron dynamics in both quantum wires— and 
quantum-dot superlattices— , where the full description on the derivation of the nonlinear 
reduced Boltzmann transport equations can be found. In the matrix Eq. ( |T0|) . we introduced 
the matrix elements (i'jj,{t) = ajj'{t) — aj^uif) via its components 



1-S 



'j, iN+l)/2 



^'^ 2h5k 



(5, 



ij'-i - 'Jj,j'+ 



2r, 
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j+j',N+l 



1 - Sj^ (jV+l)/2 

2r,- 



(11) 



While the net elastic scattering rate is given by Eq. (|H]), the total inelastic rate assumes the 
form Wj = I/t/"^ = Yl ^jj' with the scattering matrix given by 



J, J 



2n^clpLWe\^{\kji -kj\) 



(12) 



[l + cos(0j/ -0j)] • (13) 



We have /" = ff\ // = 1 - /f, n,,. = N^{\ey - e,\/n), N,{u,) = [expihcuJk^T) - 1]-^ 
is the Bose function for thermal equilibrium phonons. Furthermore, the dynamical phonon 
scattering rate Wj(t) for nonlinear phonon scattering, which is also responsible for the 
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nonlinear electron transport and electron heating due to its dependence on gj/{t), is in the 
form of 



(14) 



Once the non-equilibrium part, g'j{t), of the total electron distribution function has been 
solved from Eq. (fTOl) . the transient drift velocity, fc(i), of the system can be calculated 
according to 



It is noted that the thermal-equilibrium part of the electron distribution does not contribute 
to the drift velocity. The steady-state drift velocity of electrons is given by Vc{t) in 
the limit t — > oo. Corresponding steady-state conduction current is given by / = erii/jt'd. 
The differential electron mobility for nonlinear transport is generalized to = dv^/dJ^Q. 
Numerical simulation of these quantities in specific ANRs is the subject of the next section. 

III. NUMERICAL RESULTS AND DISCUSSION 

Figure presents our calculated electron mobilities function of applied electric 

field J-Q at T = 10 K (blue curve) and T = 6K (red curve), respectively. Clearly, from 
Fig.[2](a), a strong J^o-dependence for fie appears at a lower value of J-q at higher tempera- 
ture T in the nanoribbon. This J-Q-dependent electron mobility /Xg has its physical origins 
in the dynamical electron-phonon scattering rate (t) through its dependence on elec- 
tron the distribution function given in Eq. (1141) . Consequently, it is reasonable to expect a 
lower threshold field, J-"*, for entering into a nonlinear transport regime (J-" > J-"*) due to 
enhanced nonlinear phonon scattering at T = 10 K. The value of J-"* strongly depends on 
the parameters of the system, such as T, niD, 7o and Aq, and an analytic expression for J^* 
cannot be obtained in the nonlinear-transport regime. On the other hand, as J-q 0, fig 
is larger at T = 10 K than T = 6 K due to thermal population of high-energy states with a 
large electron group velocity. The initial decrease of fie with J-q is attributed to the gradual 




metallic 



semiconducting 



(15) 
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increase of the frictional force from phonon scattering by J-q. At T = 10 K, fie is roughly 
independent of J-q below 0.75kV/cm (linear regime). However, /ie increases significantly 
with J-'o above 0.75kV/cm (nonlinear regime). Eventually, /ie decreases with J-j) once it 
exceeds 1.5kV/cm (heating regime), leading to a saturation of the electron drift velocity. 
The electron group velocity \vj\ increases with the wave number for small \kj\ values, as can 
be seen from Eq. ([1]). However, the increase of \vj\ slows down toward its upper limit 
provided \kj\ ^ n/SW but still within the single-subband regime. The increase of /ie with 
J-Q in the nonlinear regime comes from the initially-heated electrons in high energy states 
with a larger group velocity, while the successive decrease of /ie in the heating regime comes 
from the combination of the upper limit Vj < up and the dramatically increased phonon 
scattering. In Fig.[2]^b), the calculated electron drift velocities are plotted as a functions 
of temperature when J-q = 2kV/cm (blue curve) and J-q = IkV/cm (red curve). The fact 
that /ie increases with T monotonically in both cases implies the electron scattering in two 
samples is not dominated by phonons but by impurities and line-edge roughness. Different 
behaviors in the increase of /ie with temperature can be seen from Fig.|2](b) for linear and 
nonlinear electron transport. At J-q = 2kV/cm for the high-field nonlinear transport, Vd 
(or /ie) increases with T sub-linearly. On the other hand, rises super-linearly with T 
for the low-field linear transport at J-q = IkV/cm. These different T dependence in /ie for 
linear and nonlinear transports can be directly related to the non-equilibrium part of the 
electron distribution function gj presented in Figs.|2]^c) and|2]^d). The calculated gj (left 
scaled solid curves) at T = 10 K, as well as the total electron distribution function fj (right 
scaled dashed curves), are shown in Fig.[2|^c) as functions of electron wave number kj along 
the ribbon when J-'q = 2kV/cm (blue curves) and J-q = IkV/cm (red curves). The electron 
heating may be visualized from Fig.[2t^c) when J-q = 2kV/cm as a result of thermally driven 
electrons from low to high-energy states by heat generated from a frictional force— through 
nonlinear phonon scattering. On the other hand, when J-q = IkV/cm, electrons are only 
swept by elastic scattering from the right Fermi edge to the left Fermi edge, corresponding to 
linear transport. Figure [2]^d) also allows a comparison between the calculated gj (left-hand 
scaled solid curves) and fj (right scaled dashed curves) with J^q = 2kV/cm as functions of kj 
at T = lOK (blue curves) and T = 6K (red curves). One may conclude from Fig.[2]^d) that 
the nonlinear phonon scattering becomes important at T = 10 K, while the elastic scattering 
of electrons dominates at T = 6K, similar to the explanation given for Fig.[2](c). 
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We present, in Fig.|3]^a), our numerical results for function of J-q at T = lOK 

for different electron densities riiD = 1.0 x 10^ cm~^ (blue curve) and riio = 0.8 x lO^cm"^ 
(red curve). As J-q — )■ 0, /ie is slightly increased for a higher value of tiid due to introducing 
high-energy occupied states with a large group velocity. It is further found from Fig.|ni^a) 
that a low electron density produces a small J-"* in the system. This is due to the fact that 
a small minimum energy huq ~ 2hcskp is required for initializing phonon scattering in a low 
density sample, which enhances the nonlinear phonon scattering at a chosen temperature. 
However, the enhancement of /ie due to J-q is large in a high density sample due to increased 
initial electron heating. In Fig.[3I^b), is presented at J-'o = 2kV/cm as a function of T for 
Hid = 1.0 X 10^ cm~^ (blue curve) and = 0.8 x 10^ cm^^ (red curve). Although we still 
see He increasing sub-linearly with T, as explained in Fig.[2|^b), the rate of fie enhancement 
is large for a high-rii^) sample when T is not low, where the electron heating is significant. 
This observation is complemented by the results shown in Fig.|3]^c), where electron heating 
becomes more significant for a high density sample although the nonlinear phonon scattering 
is seen to be important in both low and high density samples. This enhancement of He with 
T at high electron densities is due to more low-energy electrons being available for thermal 
excitation to occupy high-energy states through electron heating. 

The effects due to impurity scattering are compared in Figs.|l](a)-(c). In Fig.HJ^a), we 
present a comparison of mobilities as a function of J-q at T = lOK for 70 = 1.0 x 10^^ s~^ 
(blue curve) and 70 = 1.0 x 10^^ s^^ (red curve). As J-q — )■ 0, fie is greatly reduced by 
strong impurity scattering with 70 = 1.0 x 10^"^ s^^ in the linear regime. In addition, for 
7o = 1.0 X 10^^ s~^, J-"* is pushed upward from about 0.75kV/cm to a higher value at 
1.75kV/cm, leaving us with a roughly J-'o-independent fig in this case for the whole field 
range shown in this figure. The comparison for T-dependence of fie is presented in Fig.Hl^b) 
for J-Q = 2kV/cm, where a sub-linear increase of fie with T for weak impurity scattering 
is switched to a super-linear relation in the strong impurity scattering case. The effect of 
impurity scattering can also be seen from the calculated Qk and fk as functions of k at 
T = 10 K and J-q = 2kV/cm in Fig.Hl^c), where the electron heating with 70 = 1.0 x 10^'^ s~^ 
through nonlinear phonon scattering is completely masked by the strong impurity scattering 
at 7o = 1.0 X 10^^ s~^, leading to relative cooling of conduction electrons in a nanoribbon in 
the linear regime. 
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At relatively low temperatures T < 10 K and low concentration of electrons (riiD = 
lO^cm"^) in the lowest conduction band, as well as medium impurity concentration 
= 10^°cm~^, the single conduction band approximation can be justified for graphene 
nanoribbons. Indeed, the energy separation between the first and second conduction sub- 
band is Tchup/SW ~ 0.136 eV for W = 50 A. On the other hand, elastic and inelastic decay 
rate ^71 < hjo = 0.065 eV, as well as the Fermi energy eFij^nio /2) ^ 0.01 eV, are much 
smaller.— More generally speaking, the model is valid provided that the scattering parame- 
ters and Fermi energy are of the order of the subband gap Eg = 1.38/1^ eV with W expressed 
in nanometers. 

The increase of ribbon's width W will enlarge the electron group velocity Vj < vm ~ [1 + 
(2/3ni£)14^)^]~^/^, as shown by Eq. ([T]), and significantly decrease the effect of line-edge rough- 
ness scattering at the same time, as can be seen from the factor 71 ~ {j^pp I'iWY {5h/WY 
in Eq. ([9]). Fixing Aq = 200 A and making J^q — )■ 0, we do see from Fig.[5](a) that /ie 
is dramatically increased when W is changed from 40 A (71 = 1.7 x 10^'^ s^^) to 50 A 
(71 = 4.28 X 10^'^ s~^) due to enlarged group velocity of electrons and reduced line-edge 
roughness scattering. In addition, the magnitude in the enhancement of with J-q is large 
in the sample with W = hO K because the weakened elastic scattering, and therefore, the 
relatively strong initial heating of electrons by phonon scattering to high energy states with a 
larger group velocity. This is also reflected in the temperature dependence of presented in 
Fig.[5](b). However, the rate for /ig enhancement with T is almost the same in the two cases, 
implying a similar magnitude of nonlinear phonon scattering of electrons in both samples. 
This can be very well justifled by the calculated gj shown in Fig.[5](c), where the increase of 
W from 40 A to 50 A only slightly enhances the asymmetry of gj with respect to ±fcj, but 
the unique feature of electron heating is largely unchanged. 

Finally, the effect of correlation length for the line-edge roughness on the transport is 
demonstrated in Figs.[6](a)-(c) by flxing = 50A and changing Aq from 200 A to 50 A. As 
shown in Eq. ([9]), the line-edge roughness scattering can be either reduced or enhanced by 
decreasing Aq, depending on \kj\ <^ 1/2 Aq or \kj\ ^ l/2Ao. For our chosen sample with 
riiD = 1.0 X lO^cm"^, we flnd that the condition \kj\ ^ l/2Ao is satisfled in the low-fleld 
limit {\kj\ ~ k-p), while \kj\ ^ l/2Ao holds for the high fleld limit due to electron heating. 
Therefore, we flnd from Fig.[6t^a) that /ig is increased as J-q — )■ when Aq is reduced to 50 A in 
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the low field regime. However, the value of T* for is pushed upward for Aq = 50 A in the 
high-field regime, which is further accompanied by a reduced magnitude in the enhancement 
of [If. with Tq. This anomalous feature associated with reducing Aq also has a profound 
impact on the T-dependence of /ie as shown in Fig.[6]^b), where the increasing rate of /ie 
with T in the high field regime becomes much lower with Aq = 50A than for Aq = 200 A. 
In addition, the calculated gj in Fig.[6](c) displays an anomalous cooling behavior (with a 
smaller spreading of jj in the kj space) in the high field regime for conduction electrons as 
Aq is reduced to 50 A, similar to the results shown in Fig.Hl^c). 

We have performed similar calculations with metallic^^ ribbons of comparable width 
= 50 ± 5 A and found no perceptible difference in the nonlinear electron transport and 

heating provided that the rest of the parameters (electron and impurities concentrations and 

decay rates) are kept the same. 

IV. CONCLUDING REMARKS 

In conclusion, we have found there is a minimum electron mobility just before a threshold 
for an applied electric field when entering into the nonlinear transport regime, which is 
attributed to the gradual build-up of a frictional force from phonon scattering by the applied 
field. We also predict a field-induced mobility enhancement right after this threshold value, 
which is regarded as a consequence of initially-heated electrons in high energy states with 
a larger group velocity in an elastic-scattering dominated graphene nanoribbon. Moreover, 
we have discovered that this mobility enhancement reaches a maximum in the nonlinear 
transport regime as a combined result of an upper limit for the group velocity in a nanoribbon 
and a field-induced dramatically-increased phonon scattering in the system. Finally, the 
threshold field can be pushed upward and the magnitude in the mobility enhancement can 
be reduced simultaneously by a small correlation length for the line-edge roughness in the 
high-field limit due to the occupation of high-energy states by field-induced electron heating. 

When the electron density is increased in a graphene nanoribbon, multi-subband trans- 
port occurs, the field-induced mobility enhancement is expected to be reduced, and the 
effect of electron-electron scattering needs to be included. When the lattice temperature 
becomes high, on the other hand, both the optical phonon and inter- valley scattering should 
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be considered. As the width of a nanoribbon is modulated, a periodic potential along the 
ribbon will form, leading to a graphene nanoribbon superlattice with additional miniband 
gap opening at the Brillouin zone boundaries. The results of the current research is expected 
to be very useful for the understanding and design of high-power and high-speed graphene 
nanoribbon emitters and detectors in the terahertz frequency range. 

Acknowledgments 

This research was supported by contract # FA 9453-07-C-0207 of AFRL. DH would like 
to thank the Air Force Office of Scientific Research (AFOSR) for its support. 



1 A. K. Geim and K. S. Novosolev, Nat. Mater. 6, 183 (2007). 

^ K. S. Novosolev, A. K. Geim, S. V. Morozov, D. Jiang, M. 1. Katsnelson, I. V. Grigorieva, S. 

V. Dubonos, and A. A. Firsov, Nat. (London) 438, 197 (2005). 
3 C. Berger, Z. Song, X. Li, X. Wu, N. Brown, C. Naud, D. Mayou, T. Li, J. Hass, A. N. 

Marchenkov, E. H. Conrad, P. N. First, and W. A. de Hear, Sci. 312, 1191 (2006). 
^ A. H. C. Neto, F. Guinea, N. M. R. Peres, K. S. Nonoselov, and A. K. Geim, Rev. Mod. Phys. 

81, 109 (2009). 

^ P. Avouris, Z. Chen, and V. Perebeinos, Nat. Nanotechnol. 2, 605 (2007). 
6 T. Ando, J. Phys. Soc. Jpn. 75, 074716 (2006). 

J. H. Chen, C. Jang, S. Adam, M. S. Fuhrer, E. D. Williams, and M. Ishigami, Nat. Phys. 4, 

377 (2008). 

8 H. M. Dong, W. Xu, Z. Zeng, T. C. Lu, and F. M. Peeters, Phys. Rev. B 77, 235402 (2008). 

9 N. M. R. Peres, J. M. B. Lopes dos Santos, and T. Stauber, Phys. Rev. B 76, 073412 (2007); 
also see T. Stauber, N. M. R. Peres, and F. Guinea, ibid. 76, 205423 (2007). 

° V. V. Cheianov and V. I. Falko, Phys. Rev. Lett. 97, 226801 (2006). 
^ W. Xu, F. M. Peeters, and T. C. Lu, Phys. Rev. B 79, 073403 (2009). 

2 S. Y. Liu, X. L. Lei, and N. J. M. Horing, J. Appl. Phys. 104, 043705 (2008). 

^ Y.-M. Lin, C. Dimitrakopoulos, K. A. Jenkins, D. B. Farmer, H.-Y. Chiu, A. Grill, and P. 
Avouris, Sci. 327, 662 (2010). 



13 



1^ T. Mueller, F. Xia, and P. Avouris, Nat. Photonics 4, 297 (2010). 

Y.-M. Lin, V. Perebeinos, Z. Chen, and P. Avouris, Phys. Rev. B 78 161409 (2008). 

16 T. Fang, A. Konar, H. Xing, and D. Jena, Phys. Rev. B 78, 205403 (2008). 

17 X. Wang, Y. Ouyang, X. Li, H. Wang, J. Guo, and H. Dai, Phys. Rev. Lett. 100, 206803 (2008). 

18 D. H. Huang and G. Gumbs, J. Appl. Phys. 107, 103710 (2010). 

19 D. H. Huang, S. K. Lyo and G. Gumbs, Phys. Rev. B 79, 155308 (2009). 

20 S. K. Lyo and D. H. Huang, Phys. Rev. B 73, 205336 (2006). 

21 L. Brey and H. A. Fertig, Phys. Rev. B 73, 235411(2006). 

22 Such intcr-valley scattering would require momentum transfer comparable with the distance 
between K and K' points. 

23 D. H. Huang, T. Apostolova, P. M. Alsing and D. A. Cardimona, Phys. Rev. B 69, 075214 
(2004). 

2"^ T. Ando and T. Nakanishi, J. Phys. Soc. Japan, 67, 1704, (1998). 

25 K. Wakabayashi, Y. Takane, M. Yamamoto and M. Sigrist, New J. Phys. 11, 095016, (2009). 

26 L. Brey and H. A. Fertig, Phys. Rev. B 75, 125434 (2007). 

27 S. K. Lyo and D. H. Huang, Phys. Rev. B 66, 155307 (2002). 

28 A distinct line must be drawn between equilibrium distribution function in absence of 

applied electric field and stationary solution of the transport equation lirn fj{t). 

29 h'ji < hjo is vahd for Aq/PF < 6 

3° Two carbon atoms have been added/subtracted to the semiconducting ribbon across its width. 



14 



y 




FIG. 1: Schematic illustration of the model employed in our numerical calculations. The parameters 
used in this figure are justified in the text. 
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FIG. 2: (Color online) (a) Calculated electron mobilities /ie as a function of applied electric field 
J-Q at T = lOK (solid squares on blue curve) and T = 6K (solid circles on red curve); (b) Electron 
drift velocities as a function of temperature T at J-'q = 2kV/cm (blue curve) and J^q = IkV/cm 
(red curve); (c) Non-equilibrium part of (g^, left-hand scaled solid curves) and total (fk, right- 
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